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The dynamics of an exciton-LO-phonon system after an ultrafast optical excitation in an 
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cess in the absorption spectrum and time evolution are highlighted, providing a possible way of 
experimental investigation on phonon anharmonicity effects. 
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I. INTRODUCTION 

Unlike natural atoms, semiconductor quantum dots 
(QDs) always form part of a macroscopic crystal. The 
interaction with the quasi-continuum of lattice degrees 
of freedom (phonons) constitutes an inherent feature of 
these nanometer-size systems and cannot be neglected in 
any realistic modeling of QD properties, especially when 
the coherence of confined carriers is of importance. The 
understanding of the decisive role played by the carrier 
coupling to lattice modes has increased recently due to 
both experimental and theoretical study, including the ef- 
fects of the carrier-phonon interaction in a system driven 
optically on ultrafast time scales 1 * 2 * 3 *^. The phonon- 
induced decoherence seems to be crucial for any quan- 
tum information processing application and for any nan- 
otechnological device relying on quantum coherence of 
confined carriers^. 

One of the issues of interest is the evolution of the 
carrier-lattice system after ultrafast excitation of carri- 
ers. Experiments performed on bulk and quantum well 
systemsi2i2ii2iii show strongly non-Markovian behavior 
of carrier relaxation on 100 fs timescale due to pertinent 
system memory, yielding access to the quantum kinetic 
regime of the system dynamics where fundamental ideas 
of quantum mechanics may be tested. 

For systems confined in all 3 dimensions, where purely 
longitudinal optical (LO) phonon relaxation is forbid- 
den by energy conservation, coherent LO phonon dy- 
namics may still be induced by lattice relaxation af- 
ter non-adiabatic change of carrier distribution (dress- 
ing process^). In this process, an ultrafast excitation 
of confined carriers results in the lattice state which no 
longer corresponds to the energy minimum and therefore 
is unstable. Before the excess energy may be transferred 
outside the QD area (by radiating out phonons), a tran- 
sient dynamics takes place, manifested by the oscillations 
(beats, overdamped for acoustical phonons) in the optical 
response of the systemi. 

Femtosecond time-resolved experiments showing the 
joint LO-phonon evolution (beats) in such quasi-0- 
dimensional systems were performed only on semicon- 
ductor nanocrystalsiSiiiii, where confined and surface 



modes come into play, leading to different properties of 
the carrier-phonon coupling (e.g. size dependence^). 
For self-assembled semiconductor QDs, only 1 ps 
timescale dynamics was reported 1 . However, some the- 
oretical models, developed for the description of nonlin- 
ear optical response in bulk or 2D systemsiSiii, may be 
successfully applied to the confined (0D) case with dis- 
crete carrier spectrum. Within such independent boson 
models it was possible to describe the linear^i^^ and 
nonlinear— response in the limit of infinitely short pulses 
and to account quantitatively^ for the short-time behav- 
ior of optical polarization observed in the time-resolved 
spectroscopy experiment 1 . Although the observed effects 
may be attributed only to acoustic phonons 3 ", the the- 
ory predicts also LO phonon effects for sufficiently short 

Due to the weak dispersion of LO modes effectively 
coupled to the confined carriers, the dephasing of LO 
phonon beats is extremely slow (100 ps timescale)4. On 
the other hand, anharmonic damping reduces the life- 
time of a zone-center LO phonon in GaAs to 9.2 ps at 
low temperatures^. Hence, including the anharmonic- 
ity is essential for a correct description of the confined 
carrier-LO-phonon dynamics. The actual nature of an- 
harmonic phonon-phonon interactions is, unfortunately, 
still unclear. Experimental results for the bulk GaAs^S 
suggest dominant role of the decay channel to zone-edge 
LO and transversal acoustical (TA) phonons, although 
another interpretation of the same data proposes a major 
contribution of the decay into two longitudinal acoustical 
(LA) phononsSi. The latter channel has been also studied 
theoretically both for bulk 22 and confined 23 LO phonons. 
Quite unexpectedly, phonon dispersion measurements 24 
suggest that the zone-edge process may be forbidden by 
energy and momentum conservation. Since the anhar- 
monic effects seem to be essential e.g. for correct model- 
ing of relaxation phenomena— 4 2iii 2 Jii 2 i, this ambiguity is 
a real problem. 

In the present paper we describe theoretically the 
quantum evolution of LO phonons damped by anhar- 
monicity after an ultrafast carrier (exciton) excitation. 
We focus on the effects of the zone-edge decay process 
and its possible signatures in the spectral line shapes and 
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temporal dynamics. We single out various regimes both 
for single mode evolution and for the collective dynam- 
ics of all the coupled phonons. For each of these cases 
we calculate the linear optical absorption for frequencies 
corresponding to LO phonon peaks as well as the evolu- 
tion of the reduced density matrix after an ultrafast ex- 
citation. Both characteristics (in frequency and time do- 
main) are calculated within a single formalism, based on 
the Green function technique, restricted to single-phonon 
effects (which is a reasonable approximation in the weak 
coupling case at low temperatures) but incorporating 
the contribution from optically inactive exciton states 
excited by nonadiabaticity of LO phonon dynamics 2 ^. 
We show that the zone-edge decay of LO phonons mani- 
fests itself clearly in the spectral properties of the system 
for certain energy relations between the zone-center and 
zone-edge phonons. Although these energy relations are 
a material feature under fixed conditions, in view of the 
existing uncertainty concerning the phonon anharmonic- 
ity channels it seems interesting to study all possible sce- 
narios. Moreover, even if the zone-edge process turns 
out to be forbidden at normal pressure, it is likely to 
appear at high pressures. Indeed, the measured pressure 
(volume) dependence of phonon frequencies 29 shows that 
the energy of the zone-center LO phonon increases with 
growing pressure, while the total energy of the zone-edge 
LO and TA phonon pair slowly decreases. The cross- 
over pressure and the phonon replica line shapes close 
to the cross-over would yield information both on the 
phonon energy relations and on the anharmonic coupling 
strength. 

The paper is organized as follows: The sectionlTTIintro- 
duces the model and discusses the formal quantities con- 
veniently characterizing the system dynamics. Technical 
derivation is contained in the section IIIII The physical 
discussion of the result is given in the section llVl The 
section Ivl concludes the paper with final remarks. 

II. THE MODEL AND THE FORMALISM 

We consider the Hamiltonian describing a single exci- 
ton interacting with LO phonons, 

H = Hx + H p h + Hj nt . (1) 

The first part is the exciton Hamiltonian, 

Hx = e n a+a n , 

n 

where e„ is the energy of the nth exciton state and a n , a^ n 
are the corresponding annihilation and creation opera- 
tors. Our discussion will be simplified by using the spe- 
cific structure of the excitonic spectrum in self-assembled 
InAs/GaAs systems^. Due to the much larger hole effec- 
tive mass, the lowest excited states are formed to a good 
approximation by exciting the hole. Therefore, the low- 
est excited states of the confined exciton may be mod- 
eled by products of the corresponding harmonic oscillator 



wave functions (~ exp[— (ri/li) 2 /2]; ground state for the 
electron and ground or excited state for the hole), with 
Gaussian widths k given in TableQ] The resulting excited 
two-particle states are either strictly optically inactive 
due to parity or almost inactive because the ground state 
electron and excited state hole wave functions are close 
to orthogonal even in absence of strict selection rules. 
The excitation energy typical for such inactive states 
is ~ 15 meV, consistent with the thermally activated 
contribution to the central line broadening measured in 
experiment. The lowest optically active state involves 
excitation of both the electron and the hole and its en- 
ergy is typically much higher (~ 70 meV). Thus, even the 
LO replicas of the lowest excited state(s) lie still below 
the active state. Under totally non-selective excitation 
all states contribute to the dynamics and the resulting 
evolution would be both difficult to describe theoreti- 
cally and hard to interpret experimentally. Therefore, 
we assume some degree of spectral selectivity due to fi- 
nite pulse length, sufficient to avoid the excitation of the 
optically active excited state, and include in our descrip- 
tion only the lowest, optically inactive excited state. 

The phonon Hamiltonian H pn describes the LO and 
TA phonon branches coupled by anharmonicitySS, 

#ph = 51 ^fcfrfcfrfc + ^2 w fe c fe c fe 

+ Y / W(k,k')b k b{_ k ,cl+U.c., 

k,k' 

where bk,b\, and Cfc,c^. refer to LO and TA phonons 
with momentum fc, respectively, and f2fc, cok are the cor- 
responding energies. To simplify the notation, we put 
h = 1, &b = L The anharmonic term describes a decay 
of LO phonon into a TA phonon and another LO phonon, 
which may be allowed by energy and momentum conser- 
vation only when the final states lie in the vicinity of the 
L point of the Brillouin zon o 20 i 24 . The interaction term 
accounts for the Frohlich interaction between confined 
charges and the lattice polarization induced by the LO 
deformation, 

flint = — L ^2 F 7hn ,(k)a+a n , (bk + bi^j , (2) 

where the coupling constants for the confined carriers are 
given by 

«-*<*> = -li/ss / (3) 

x*;(r c ,r,,) ((""--e*") *„.(r c ,r h ). 

where r c ,r n denote the coordinates of the electron and 
hole, respectively, eo is the vacuum dielectric constant, 
and the other elements of the notation are described in 
Table [U along with values (corresponding to InAs/GaAs 
system) used in the calculations. 

The system properties may be experimentally studied 
in the frequency or time domain. Although usual spectral 
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Effective dielectric constant- 


s 


62.6 


Optical phonon energy- 




36 meV 


Unit crystal cell volume 


V 


0.044 nm 3 


r point LO dispersion parameter 24 


/^LO.r 


-0.04 meV-nm 2 


L point LO dispersion parameters 24 : 






parallel to the zone boundary 


A-o,|| 


0.17 meV-nm 2 


perpendicular to the z. b. 


0LO,± 


0.35 meV-nm 2 


L point TA dispersion parameters—: 






parallel to the zone boundary 


0TA,\\ 


0.16 meV-nm 2 


perpendicular to the z. b. 




-0.04 meV-nm 2 


Zone-center phonon lifetime- 


TLO 


9.2 ps 


Wave function widths: 






electron, in-plane 


lc,± 


4.0 nm 


hole, in-plane 


ih,± 


3.3 nm 


electron and hole, z-direction 


I, 


1 nm 



TABLE I: The parameters used in the calculations. 



functions describe the equilibrium properties of the sys- 
tem and are relevant only in the linear response regime, 
in the limiting case of ultrafast excitation they are suffi- 
cient for the description of the system evolution also in 
the strongly driven case. This is possible due to sepa- 
ration of the time scales (external driving vs. phonon 
response): the laser pulse prepares an initial state while 
the subsequent evolution consists in relaxation to equi- 
librium. 

To see this more formally, let us note that the elec- 
tromagnetic wave is coupled only to carrier degrees of 
freedom and not to lattice modes so that all information 
accessible to optical measurement is given by the reduced 
density matrix of the carrier subsystem, 



(*)= Tr L [U(t)p(0)tf(t)}, 



Q 



where p(0) is the initial state of the total system and 
the trace is taken over the lattice (phonon) degrees of 
freedom. Assuming that the state n is the only relevant 
optically active state in the frequency sector of interest, 
the optical properties are determined by the non-diagonal 
element of the reduced density matrix of the carrier sub- 
system, proportional to the optical polarization, 



g n (t) = (vac|e(t)|n) = Tr [p(0)a n (t)] 



(4) 



where |vac) stands for vacuum (empty dot). Let us con- 
sider a change in the carrier subsystem performed by 
the ultrafast pulse that takes the initial vacuum state 
instantaneously into a certain superposition of the vac- 
uum state and the single-exciton state n, i.e., \ip) = 
cos Tf\va,c) + sin ^ajjvac), where a is the area of the ul- 
trafast driving pulse. Due to its inertia, the lattice is left 
in its initial state, i.e. in thermal equilibrium described 
by the density matrix p^. The further dynamics, after 
switching off the pulse, is generated by the Hamiltonian 
(|T|). conserving the number of excitons. Inserting the ini- 



tial state p(0) — \ip){^\ ® Pl into Eq. (@J one gets 

1 

g n = - smal n (t), 

where the exciton correlation function is defined by 

J n (t) = (a n (i)4(0)) (5) 

((■) denotes the thermal average). Thus, the arbitrary 
angle a factors out and the evolution is described uni- 
versally by a single function I(t). 

The function corresponds to the "survival ampli- 
tude", i.e., to the overlap between the state of an exci- 
ton created instantaneously at t = and evolving under 
the full carrier-phonon Hamiltonian until the time t and 
the unperturbed exciton state. For our discussion it is 
essential that this correlation function is related to the 
standard Green functions formalism. Indeed, defining 
the retarded Green function 



G nn ,(oj) = J dte iut e(t)({a n (t),al}) 

and the corresponding spectral density 
A n {u) = -21mG„„(w), 

one can write^i 

/ \ -i f°° 

A n (u) = [l-e-"l T ) / I n {t)e wt dt 

J — OO 

I n {t)e lult dt. 



(6) 
(7) 



The final approximation is justified by the fact that for 
a typical semiconductor, the exciton energy is 1-2 or- 
ders of magnitude larger than the thermal energy even 
at room temperature, so the Boltzmann factor in this 
formula may be neglected. The experimentally measur- 
able absorption is proportional to A n (ui). Hence, both 
the linear spectrum and the evolution after a (possibly 
strong) ultrashort pulse is conveniently described within 
the Green function approach. 



III. EXCITON GREEN FUNCTION 

In order to find the diagonal Green function G nn (iv) we 
write the Dyson equation in the Matsubara formalism, 

>(ip„)Gi> m (ip„). 

In order to approximately invert this system of equations 
we note that the diagonal Green functions are of zeroth 
order in the carrier-phonon interaction while the non- 
diagonal ones, as well as the the mass operator elements 
S n „', are at least of the second order. Hence, the system 
may be inverted to the leading order to yield the diagonal 
exciton Green function 



G n {ip„) = G nn {ip v ) 



1 



IPv 
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where E n (ip u ) = E nn (ip v ). 

Throughout the paper, the Greek indices will be used 
for numbering Matsubara frequencies, while the Latin 
ones for exciton states (including all the relevant quan- 
tum numbers). In the single-exciton approximation, the 
exciton may be treated as a fermion, with the correspond- 
ing frequencies p v . The bosonic frequencies are denoted 
by lo^, while e n denotes the bare exciton energy in the 
eigenstate n. 

The corresponding retarded functions are 

G n (u)) — G n (ipv)\i Pv ^uj+iS 
E n (u) = E n (ip u )\i Pu ^ u+ is. 

Writing the retarded mass operator as E n (ui) — A n (ui) — 
Z7 n (w), one has 



A n {ui) + ijn{u) + iS 



(8) 



The energy levels E n for the exciton interacting with 
LO phonons (the confined excitonic polaron2i22) may be 
found as the poles of this function; their real parts are 
given by the equation 



E. 



n c n 



e„- A n (E n ) = 0. 



(9) 



Let us note that the correction A n (u>) in the denominator 
of Eq. |JSJ is important only near e„, where the other 
term vanishes (otherwise, it is a higher-order correction). 
Therefore, it may be replaced by 



A„(cj)« A n (E n ) + (u - E n ) 



dA n (u) 
du> 



and combined with e„ to give E n , in accordance with 
Eq. ©■ 

Thus, one may write 



G n (w) = Z n 



1 



U) - E n + «7n(w) ' 



(10) 



where 



Restricting the discussion to the weak coupling case, 
we are interested in the single-phonon approximation, ex- 
cept for including the lowest-order anharmonicity effects 
into the LO phonon propagation. Therefore, we neglect 
the vertex correction in the mass operator E n (ip v ) and 
replace the exact phonon Green function by the free one 
perturbed by the coupling to TA phonons in the leading 
order. The mass operator in this approximation is given 
by 



E n (ip„) = 

-l§jf E E \Fnn>(k)\ 2 G n ,{w» - i^)D(k, . 
n' k n 

and the LO phonon Green function is 



(12) 



D(k,ico fJ ,) 



-(iu^) 2 + n 2 + 2n k Tl(k, iujf,) ' 



(13) 



The polarization operator in the proposed approximation 
involves the lowest order anharmonicity correction 



u=E„ 



n(fe,KJ M ) 



44 E \W(k,kTD { Lkk-k\i^-i^)D^l(k\i^). 



(14) 



Inserting the free phonon Green functions into Eq. I|14|) and performing the summation over frequencies one gets 



n(fe, {up) 



^Ei^( fe > fc ') 



(n k > - jVfc-fc')(^lfc-fc' - ujk') {n k > +N k - k > + l)(^fc_fc/ +u>k') 
(iw M ) 2 - (O fe _ fe / - uj k <) 2 



where N k and n k are the Bose distributions for LO and TA phonons, respectively, with momentum k. 
The Green functions may be written in the representation 



D(k,iai tl ) 



dw B{k,w') 



2tt iuj^ — to 



7 , G n ,(iu) IJ ) = 



duj A n i{u>') 
lit iujfi — u' ' 



where the phonon spectral density is defined as 

B(k,oj') = -2InxD(fc,a/). 
The frequency summation in Eq. (|12J) may now be performed, leading to 



, s 1 v^, „ Nl2 f°° duj 1 f°° dw 1 ' , n ,B(k 



] B{k 1 u/)A I A^)_ 
+ u" — ip v 



(15) 



(16) 



(17) 
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From the above, using Eq. (|16fl . the retarded mass operator reads 

£„M = --]T|F lm ,(fc)| 2 / ^[n B (^) + l]S(fe,^)G„'(^-^), 

/ i J —OQ 

and, inserting Eq. (fTUjl . 

A «H = -^Y.Wnn'{k)\ 2 J ^[n B {J) + l]B{k^') 

7„m - s e i'- wi" /_ ^i»bK) + _ „, _ £J1 + 7> _ „ 



u> — u>' — E„/ 



7„(w - u/) 



(18) 

(19a) 
(19b) 



To the leading order, A n (w) and 7 n (w) may be obtained from Eq. i|19al b) by setting 7„(w) — » on the right-hand 
side, which leads to 



A„(w) 



1 /*OC j / 



B(k,uj') 

to — u' — E, n 



= 2^^|F n „'(fe)| 2 B(fe,^-i;„0[nB(^-^nO + l]- 



(20a) 
(20b) 



n' M 



According to Eqs. (|17|i and (|13fl . the phonon spectral density is given by 

8n|lmII(fc,cj) 



B{k,u) 



[to 2 - O 2 + 2H fc Rc H(fe, w)] 2 + [2f7 fc Im H(fe, u)} 2 ' 



(21) 



The essential contribution to A n (w) may be obtained by setting U(k,uj) — > in Eq. (|2*T|l : after a simple integration 
one gets 



k.n' 



l + N k 



A,, 



Hence, from Eq. I|ll|) . 



h.n' 



l + A^fe 



(£„ - En, - n k ) 2 (E n -CO- E n , + fl k ) 2 



(22) 



To find 7n(w) we get the retarded polarization operator n(fe,a>) from Eq. (|15fl by io; u — > lu + z<5; thus 

X (fc,w) = Imn(fc,w) (23) 

k' 

+(n k > + N k + 1) [<S(cj - fl k -k> - w fe /) - 5(w + Ofc-fe' + ^fe')]} ■ 

I 



According to Eq. J2U, the contribution from LO 
phonons to the spectral density B(k, Lu—E n > ) in Eq. I|20b|) 
is peaked around to — E n i w ±f2fc. Hence, the Bose func- 
tion takes the values of ne(Ofc) + 1 and iib(- ilk) + 1 = 
— rift(Q k ). If we assume that the temperature is well be- 
low 400 K, the former is close to one and the latter is 
zero. Thus, lu w E n i + Q k and we may approximate, ne- 
glecting the LO phonon energy shift due to anharmonic 



coupling, 

7«H = (24) 

CMI 2 x(k,ui-E n ,) 

N fa l * nn ' Wl [lu - E n , - n k } 2 + X 2 (k,Lo - E n .) ■ 

There are a few effects contained in this formula. First, 
it generates a broadening of the zero-phonon spectral 
peak at lu = E n due to LO phonon-assisted transi- 
tions. In the absence of anharmonic coupling (x — > 0) 
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an energy conserving Dirac S appears on the right-hand 
side, restricting the process to exact resonance. The an- 
harmonicity releases this restriction and leads to two- 
phonon anharmonicity-induced processes whenever the 
energy spacing matches an LO+L A or LO+TA pair, as 
described in detail in Ref. Il5l26l Further broadening 
comes from the usual second-order two-phonon processes 
which are beyond the leading-order approximation used 
here. For the ground state, however, all these processes 
are of phonon absorption type and are negligible at low 
temperatures. The second effect, which is of major inter- 
est to us in the present work, is related to the line shape 
of the LO phonon replica (around lo ~ E n > + Q k ) and its 
effect on the system evolution. 

In this spectral region only the first and third terms 
from Eq. (|23[1 contribute to Eq. (|24|) since the energies 
of LO phonons lie entirely above the TA branch. Out 
of these two, the first one vanishes at low temperatures. 
Moreover, since coupling to short-wavelength modes is 
favored both by the general form of the coupling (~ yk' , 
see Ref. l33j) and by the growing density of states, this 
contribution (for low fc, as selected by the exciton-LO- 
phonon form factosii^) is peaked around to — E n / « 20 
meV (for GaAs, see Ref. |2J for exact dispersion curves) 
and gives smaller contribution due to large denominator 
in Eq. Ij24(l . The remaining term corresponds to the LO 
phonon decay with emission of another LO phonon and 
TA phonon. It is large, when the total energy of the re- 
sulting phonons is close to the T-point LO phonon energy, 
i.e. near the L point in the Brillouin zon o 20 i 24 . Thus, we 
approximate the coupling constant and the Bose distri- 
bution by the corresponding L-point values, assuming 
also k near the zone center. Replacing the summation 
by integration we obtain, after these approximations, 



vir 



\W(0,k L )\ 2 (n hL +l) 



J d 3 k'5(uj - fl k -k' - Wfe')> 



where v is the unit cell volume. The problem is thus 
reduced to finding the appropriate two-phonon density 
of states. We assume that LO and TA phonon energies 
may be approximated by quadratic forms of momenta 
around each L point, given by the matrices /3lo an d /?tA' 
respectively. The index £ = 1 ... 4 numbers the four pairs 
of opposite hexagonal walls of the Brillouin zone (note 
that the four tensors differ only by orientation). Near 
the ith L point, the minimum total energy of a LO+TA 
pair with total momentum k is 



where $® = ffi 



k(3 



y LO PhA' ft ^ s c l ear that the LO phonon 
with momentum k may decay via this process only if 
> £fc = niinj^i^. One finds, 



d 3 k'S(uj - Q 



fe-fc' 



E 



dS' 



s[ 2|/3Wfc'| 



where the surface S[ in the k' space, composed by all 
states satisfying energy and momentum conservation, is 
defined by ^ + kf3^k = to (the integral is zero if this 

equation has no solutions, i.e. for lo < £i ). Thus, the 
finite bandwidth of phonon branches leads to a (lower) 
frequency cutoff in the two-phonon density of final states 
for the anharmonic process. 

For GaAs, the tensors (3^ may be reasonably approx- 
imated isotropically (and hence become equal). Defining 
(3 = (det/3^) 1 / 3 , one finds explicitly 

X (k,u>) = W*)(fc,c), 



where 



X ^{k 1 uj) = —\W{0 1 k L )\ Z (n kL 



l)/3 



-3/2 



for lo > £u and x^\k,uj) — otherwise. Since, most 
probably, the zone-edge LO+TA process is only one of 
many LO phonon damping mechanisms, we shall phe- 
nomenologically add to x(k,uj) an additional contribu- 
tion xoj corresponding to all the other damping pro- 
cesses. 

As discussed above, 7 n (w) contains various contribu- 
tions, related to the broadened zero-phonon line, LO 
phonon replica, and acoustic phonon sidebands^ (not in- 
cluded in our discussion). As long as these features are 
spectrally separated, as is the case for the ground exci- 
ton state at low temperatures, they manifest themselves 
by processes on different timescales ( 1 ns exponential 
decaji^, 1 ps overdamped acoustic phonon dynamics 2 , 
100 fs LO phonon beats damped on 10 ps timescale as 
discussed below). These processes may, therefore, be ex- 
perimentally separated even if they are superposed on 
one another. As we are mostly interested in the LO 
phonon features we simply model the low-frequency sec- 
tor by a Lorentzian, allowing for a certain broadening of 
the zero-phonon line. Separating this low-frequency part 
from the peak around lo = Eq + fl k , we write the spectral 
density for the ground exciton state in the form 



A (lo) 



2Z n 



lo 



(Lo-E r+^ 



(25) 



77 l2 $0n'( fe )9fc(^ - fife - E n '), 



n'k 



where 



$nn'(*0 



5fc(w) 



|-Fnn'(fc)| S 



{E n i — E n + flk) 2 



x(k,Lo + Q k ) 
+ X 2 {k,Lo + il k )' 



(26) 



(27) 



and 70 accounts for the broadening of the zero-phonon 
line. The first contribution corresponds to the central 
peak while the second part describes the exciton dressing 
with a continuum of LO phonon modes. 
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IV. THE TIME EVOLUTION: DAMPED LO 
PHONON BEATS 

The time-dependent correlation function for t > 0, de- 
rived from Eq. i|25|) according to Eq. J7J, is 

I (t) = (28) 

n' ,fc 

i 



where 

9h{t) = J ^g k (u)e-^ = \g k (t)\e iVh{t) . (29) 

Using Eq. 12L'I) in the low temperature limit (N k — + 0), 

one gets 



IM*)I* = 1- 



2^ *°"( fe ) I + 2 a? Z $ 0n(fe)|5fe(*)l cos [(£J - E n - Q k )t + <p k (t)] . 

n,k J n.k 



(30) 



Hence, the evolution of each LO phonon mode in the 
polarization cloud may be treated separately. The total 
effect emerges as a superposition of such individual con- 
tributions, weighted by the coupling strengths. It should 
be noted that in spite of the growing energy separation 
in the denominator of the coupling strength l|2t)[l (for 
n = 0), the coupling to higher levels (n' > 1) is stronger 
than that to the n' = level because of the charge can- 
cellation which reduces the diagonal coupling constants 
Foo{k)k 

Let us study a few examples of the system evolution. 
The numerical results presented below have been ob- 
tained with the GaAs parameters shown in Tab. |U The 
anharmonic coupling strength | W(0, fci,)| was determined 
from the requirement that the zone-edge LO+TA decay 
process should be responsible for 90% of the experimen- 
tally known decay probability of the T-point LO phonon. 
A specific LO mode fc is damped by the zone-edge pro- 
cess if x(fe, Ofe) > 0, i.e. if its energy lies over the density 
of states (DOS) cutoff, il k > £ k . In our parabolic disper- 
sion model one has 




35.5 36 36.5 200 400 600 800 
hco [meV] t [ps] 



Q k — £fc — 

^0 — ^fcr 



U k r 



max, k 



FIG. 1: Single mode line shape (left) and temporal evolution 
(envelope of LO oscillations, right) for three modes with wave 
vectors as shown, for tto — Qk, — ^k, = —0.5 meV. 



where the last term is negative except for a very small 
positive contribution (sensitive to the crudeness of disper- 
sion modeling) in a narrow solid angle around the (111) 
and equivalent directions (towards the L points). Hence, 
if the decay process is forbidden for the T-point mode it is 
also forbidden for a vast majority of all the modes. Two 
different regimes may arise when this process is allowed 
for a given mode: 

If the decay results in states distant from the zone edge 
then the density of final states is also far from its edge. 
Quantitatively, for 



V^k - Ik > 



v\W(0,k L )\< 
2tt/3 3 / 2 



("fe L + !)> 



one has 



duj 1 



«1, 



and the value of x(fe, is almost constant around the LO 
phonon energy, leading to the Lorentzian line shape (for 
this specific mode) and exponential decay (FigOJi). The 
decay time of oscillations induced by a mode in |/„(£)| 2 is 
twice longer than the decay time of the mode occupation 
(which would be governed by |fffc| 2 ). In fact, the homo- 
geneity of the "regular" decay (far from the DOS edge) 
justifies representing all such decay processes by a single 
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constant xo- On the other hand, if the anharmonic decay 
process leads to zone-edge states, then £ fc is close to £l k 
and the line shape is determined mostly by the frequency 
dependence of w ) and becomes non-Lorentzian (Fig. 

If the mode energy lies below the DOS cutoff, the mode 
is damped only by the other processes, represented by xo 
(Fig. (IJ;). The vicinity of the DOS edge still results in 
dephasing effects (by virtual processes) manifested by a 
non-Lorentzian line shape. 

In principle, all the three scenarios may be present 
among modes forming the polaronic cloud, leading to 
intricate system dynamics. However, since it is rather 
impossible to excite a single, selected mode, the observ- 
able evolution always results form superposing the dy- 
namics of all the modes forming the polaron cloud. Such 
dynamics (and the corresponding line shape) averages 
the contributions from different modes. Moreover, the 
joint evolution is affected by the dephasing effect result- 
ing from dispersionSdiiS. Let us note that, according to 
Eqs. I|25I28[) . several contributions may be present in the 
system evolution, related not only to the ground state 
itself but also from non-adiabatic excitation of optically 
inactive excited states, the latter being much more pro- 
nounced due to lack of charge cancellation. However, 
the excited state contributions correspond to higher fre- 
quencies and may be selectively turned off by modifying 
the pulse duration (i.e. its spectral width). Below we 
present the results assuming the excitation of the ground 
state only and of the ground and first excited state. We 
neglect the broadening of the central line which, in the 
time range of interest, would produce an additional weak 
slope in the time evolution. 

Before discussing the results, let us note that in the 
case of a single exciton state and without any phonon 
anharmonicity the problem is reduced to the independent 
boson model which is solvable exactly 2 . This allows us 
to estimate the accuracy of our approximate approach by 
comparing it with the exact result in this limiting case. 
As shown in Fig. [5^ the Green function technique gives 
very exact results, which is due to the fact that the overall 
phonon-induced perturbation is very weak. 

The results of our calculations for the full case includ- 
ing anharmonic phonon decay processes are shown in 
Fig. [2>d. We calculated the shapes of the LO replica 
lines for various possible energy relations at the zone 
edge. If the momentum and energy conservation for all 
effectively coupled modes (with spectral width Att^) is 
satisfied by final states lying far from the zone edge, i.e. 
if f2o — Ofc L — Wfc L 3> Aflk (large energy overlap at the 
zone edge) , the difference of decay times between modes 
becomes small and beats corresponding to all the modes 
decay with the same lifetime r, equal to the doubled life- 
time of the bulk zone-center LO phonon tlo ■ In this case, 
the damping dominates over dephasing. However, since 
the damping and dephasing times are still comparable, 
the dispersion-induced dephasing effect is not negligible 
and reduces the decay rate of the whole phonon dressing 




h(0 [meV] t [ps] 

FIG. 2: (a) The hypothetical line shape of the LO phonon 
replica of the ground exciton state without anharmonic damp- 
ing (broadened only due to dispersion) and the corresponding 
envelope of LO phonon beats (line) compared to the exact re- 
sult (squares), (b-d) The line shapes of the phonon replicas 
corresponding to the ground state (dashed lines) and to the 
first excited optically inactive state (solid lines) and the corre- 
sponding envelope of the phonon beats for three values of the 
band overlap as shown. The central frequency luq corresponds 
to E n + Qq for each line. In plots (b) and (c) it is assumed 
that the zone-edge process contributes 90% of the decay prob- 
ability. In (d), where the zone-edge process is forbidden, the 
anharmonic LO-TA coupling is the same as in (c). 



cloud oscillations below the individual mode rate 2tlo 

(Fig. Eh). 

If, on the other hand, the zone-edge energy overlap is 
small, then some modes will be weakly affected by the 
zone-edge mechanism and some others will be unaffected 
by it. Thus, some modes off the T point live much longer 
than the zone-center mode and the overall decay time 
of the phonon beats may become longer than 2tlo (Fig. 

Et). 

Finally, let us consider a situation when the LO+TA 
process is forbidden for the zone center phonon. Al- 
though some other modes still may decay via this pro- 
cess, this effect is too weak to be noticeable in the overall 
dynamics and the phonon replica line again takes the 
Lorentzian shape with decay time shorter than the tlo 
due to dispersion. 
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V. CONCLUSIONS 

In this paper we have studied the quantum evolution of 
the exciton-LO-phonon system in an InAs/GaAs quan- 
tum dot, focusing on the possible signature of the zone- 
edge LO^LO+TA anharmonic decay in the line shapes 
of the LO phonon replicas and in the time evolution of 
the system. Apart from including anharmonic damping, 
we have further extended the description by incorporat- 
ing contributions from optically inactive excited states. 
We have shown that the vicinity of the edge of density of 
final states leads to the appearance of long-living modes 
and to a deviation from the symmetric Lorentzian line 
shape. In such case, the decay time of the LO phonon 
beats becomes longer than 2tlo> where tlo is the mea- 
sured decay time of the zone-center LO phonon popu- 
lation. This situation appears when the total energy of 
the L-point phonons (LO+TA) is only slightly lower (by 
~ 0.1 meV) than the T-point LO phonon energy. On the 
other hand, when this energy overlap becomes larger, the 
line shape becomes very close to Lorentzian and the de- 
cay time of the LO phonon beats is shorter than 2tlOj 
due to additional contribution from dispersion-induced 
dephasing. 

In view of long coherence time of the undamped LO 
phonon oscillations (due to weak dispersion), the line 
shapes of the phonon replicas and the time evolution 
of beats is always determined by the anharmonic damp- 
ing processes. Therefore, the damping effect must be 
included in any description of carrier-LO-phonon dy- 
namics. Moreover, in some special cases, the detailed 
structure of phonon bands may noticeably influence the 
measurable properties. It may therefore be possible to 
extract information about the anharmonic LO phonon 
damping channels from the QD spectroscopy data. Al- 



though the phonon dispersion is obviously a material 
characteristics, it may be modified by changing the mate- 
rial composition and experimental conditions (e.g. pres- 
sure) , so a range of dynamical regimes might be observed. 

It may be expected that an actual experiment will be 
performed on the ensemble of dots by nonlinear optical 
techniques, e.g. four- wave- mixing (FWM). It is known 
that the time-resolved FWM results cannot be repro- 
duced exactly by a linear spectrum or by a single pulse 
dynamics 3 . Unlike for the Markovian dephasing effects, 
in the case of the lattice response related to pure dephas- 
ing after a fast excitation there is no simple relation be- 
tween the spectra and time-resolved evolutions obtained 
in linear and nonlinear experiments. However much simi- 
larity still remains, due to the unique underlying physical 
effects. For instance, the variation of damping times and 
the non-Lorentzian shape of the LO replica reflect the 
structure of the allowed final states for an anharmonic 
decay process and may be expected to appear in any 
kind of experiment. Hence, the results obtained here are 
at least qualitatively relevant also for a nonlinear experi- 
ment, although an exact extension in this direction would 
certainly be desirable. Another important generalization 
should consist in including final-length pulses^* 3 ^^. 
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